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Abstract 

We introduce a deformation of the method of characteristics valid for Hamil- 
tonian perturbations of a scalar conservation law in the small dispersion limit. 
Our method of analysis is based on the 'variational string equation', a functional- 
differential relation originally introduced by Dubrovin in a particular case, of which 
we lay the mathematical foundation. 
(-h ' Starting from first principles, we construct the string equation explicitly up to the 

O ^ fourth order in perturbation theory, and we show that the solution to the Cauchy 

problem of the Hamiltonian PDE satisfies the appropriate string equation in the 
small dispersion limit. We apply our construction to explicitly compute the first 
two perturbative corrections of the solution to the general Hamiltonian PDE. In the 
KdV case, we prove the existence of a quasi-triviality transformation at any order 
and for arbitrary initial data. 

Introduction 

The present paper is devoted to the study of Hamiltonian perturbations of the scalar 
quasihnear conservation law 

u t = a(u)u x , (1) 

where u = u(x,t), i,( 6 R, and a is a non-constant function of u. This equation 
admits a formal Hamiltonian representation, given by 

Ut (x) = {u(x),H%}, H%:= I ' h(u)dx, h" = a, 

where the Poisson bracket is the Gardner-Zakharov-Faddeev bracket: 

s „ rn [SHdSK 

{H > K}:= h^T x ^i dx ' (2) 
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for every pair of functionals H, K . The variational derivative 4- is defined for any 
local functional H = J h(u, u Xl u XXl . . . )dx through the Euler-Lagrange operator £ , 
namely: 

k>0 

The standard method of characteristics provides local classical solutions of the 
Cauchy problem for ([1) by solving the implicit equation 

x + ta(u) = g'(u), (4) 

where g' is a local inverse of the initial data. The time lifespan of the solutions for 
generic initial data is finite: there exists a i c < oo, known as the critical time, be- 
yond which no classical solution exists. The point (x c ,t c ) where the solution breaks 
is known as the critical point of the solution. 

The class of Hamiltonian perturbations of equation (TTJ) considered in the present 
work is obtained by keeping fixed the Poisson bracket @, and formally deforming 
the Hamiltonian H? in the following way: 

H h := Hi + e 2 H 2 + e 4 H A + . . . . 

Here, e is a (small) parameter, every functional H k is required to be local, and the 
corresponding density to be a homogeneous polynomial in the derivatives j = 
1,2, ... , of total degree k, where degM J = j. As an example, the KdV equation 
Ut = uu x + s 2 u xxx fits in this class by choosing h = H 2 = — ^ J u 2 dx, and 
H k = 0, k > 2. 

The above class of Hamiltonian perturbations has been studied in detail by 
Dubrovin in [3], where in particular he provides a complete classification (up to 
order e 4 ) of families of commuting Hamiltonians, and he formulates an important 
universality conjecture on the critical behaviour of Hamiltonian perturbations. Ac- 
cording to the conjecture, for any perturbed PDE and for generic initial data the 
solution of the Cauchy problem exists up to some time, bigger than the critical 
time of the unperturbed equation. Moreover, around the critical point, every such 
solution is locally described by a particular solution of a 4 t,l -order ODE of Painleve 
type. The universality conjecture was proved to hold for any equation of the KdV 
hierarchy, for particular classes of solutions (5J [3]. One of the key steps for the 
derivation of the universality conjecture obtained in [4] is the introduction, for a 
special subclass of perturbations, of a formal identity, which in the KdV case takes 
the form 

x +ut=g' + e 2 {^9 {i) u 2 x + g"'u x )j + e* (^-g^u 2 x u xx 

+^9 {5) uL + !</ ""<<, <>,,, + \g [i) u xxxx - ^g {7) <) + 0(s 6 ) . (5) 

This equation clearly reminds of the characteristics equation ([1]): the aim of the 
present paper is to construct an analogue of the above formal identity for a more 
general class of Hamiltonian perturbations, and to rigorously prove its range of va- 
lidity. It turns out that both problems are non-trivial. As an application, we use 
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the new equation to explicitly find the first two corrections of the semiclassical ex- 
pansion of the solution to the perturbed PDE. It should be noted that in the case 
g = u n , n € N, then identity reduces, modulo higher order terms, to the string 
equation of the KdV hierarchy [151 113| . By analogy, and considering that our pro- 
cedure to obtain the equation follows a variational formulation, the new equation to 
be constructed will be named 'variational string equation'. 

The paper is divided into two main parts, one more algebraic in nature, and 
a second one containing analytic results. After reviewing some known facts on 
Hamiltonian perturbations in Section [TJ we consider in Section [5] the construction 
of the variational string equation. This is done according to the following scheme: 
first, we require the string equation to describe the stationary points of a certain 
(approximate) conserved quantity, namely to be of the form 

fi = o, 

OU 

for some suitable functional Sf. Besides this, at e = we want to recover the char- 
acteristics equation and at t = the equation has to provide the desired initial 
data. While the latter two conditions are very natural, the motivation for the first 
one comes from the following result, obtained by Lax in [9pl : given a (nonlinear) 
evolutionary equation, together with a conserved quantity, then the space of sta- 
tionary points of the conserved quantity is invariant under the flow. As shown in 
Section [TJ equation Q fits in this picture; we require the variational string equation 
to be of the same nature. 

As a byproduct of the above construction, we obtain a characterization of the 
Hamiltonian perturbations admitting a variational string equation: at order 0(e 2 ) 
the string equation exists for every perturbation, while at order 0(e 4 ) we find that 
a certain functional constraint has to be satisfied. On the other hand, if for a given 
perturbation a string equation exists (at some order), then it is unique. 

In Section [3] we describe how to use the variational string equation in order 
to provide approximate solutions - for small times - to the perturbed equations. 
Given the solution u of a Cauchy problem for an Hamiltonian perturbation, with 
e-independent initial data, then at a formal level the coefficients v z ,i > 1, of the 
semiclassical expansion 

u = v° + e 2 v 1 + e i v 2 + ... 

can be written in terms of the solution of the unperturbed equation v° (with same 
initial data). If the Hamiltonian perturbation admits a string equation at order 2fc, 
then the first k coefficients can be explicitly found by elementary manipulations. 
This procedure was already used in our previous paper |12j . where we wrote the 
result for the first coefficient v 1 without explaining the derivation. In the present 
work we show the details of the proof, and we compute the second term v 2 . 

Finally, Sectionals devoted to the rigorous proof of the validity of the variational 
string equation. We show that, given the solution u(x,t,e) of a Cauchy problem for 
a Hamiltonian perturbation which admits a variational string equation (at some 

An analogous result was proved independently by Novikov |14| . by considering symmetries 
rather than conserved quantities. In the Hamiltonian case the two approaches are equivalent. 
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order in e), then there exists an integer K and a domain in the (:c,i)-plane where 
the variational string equation is satisfied up to order e K , that is: 




for x and t belonging to the domain. The number K depends both on the regularity 
of the solution of the Cauchy problem and on the order in e of the string equation. 
The problem is non standard, for the functional Sf is in general ill-defined and the 
variational string equation is defined only locally. We approach it by first showing 
that the variational derivative of a formal functional satisfies a linear PDE with 
forcing term (this result, already interesting by itself, is a generalization of a result 
by Lax), and then describing in detail the domain of the variational string equation. 
After that, we are able to control the L 2 -norm of -r^-, uniformly in x and t, and to 
show that it is small in e as required. Since the value of at the boundary of the 
domain is not known, an important ingredient of the proof is the introduction of 
a family of curves in the (x, i)-plane, suitable deformations of the characteristic lines. 

The results of Section 0] are valid provided the solution of the Cauchy problem 
is regular enough in the semiclassical limit. Therefore, using the results we proved 
in |12| . they can be applied to Hamiltonian perturbations of generalized KdV type. 
In Section [5] we briefly review these results, and we prove that all the perturbative 
corrections of the Cauchy problem of KdV are rational functions of the solution of 
the unperturbed equation and its derivatives. This property, called quasi-triviality, 
was proved so far just for solutions of the unperturbed equation with never vanishing 
first derivative and for initial data with a depending on e in a prescribed way [Hill). 

1 Hamiltonian perturbation of quasilinear PDEs 

Here we briefly review the Dubrovin-Zhang construction [3] of Hamiltonian per- 
turbations for quasilinear conservation laws. We focus on the scalar case, which is 
treated in detail in @]. In order to fix the notation, we recall first some well known 
fact about the unperturbed case; we then consider the theory of Hamiltonian per- 
turbations, and finally we discuss the problem of finding approximate solutions (for 
small times) to the perturbed equations. 

1.1 Quasilinear conservation laws 

Consider the Cauchy problem for a scalar conservation law 



with a(v) any non-constant smooth function. We look for classical solutions of ©, 
that is for a bounded differentiable function, with bounded space derivative, such 
that its partial derivatives satisfy © identically on K x [0, t*[ for some positive 
time t* . The standard way to construct solutions of the Cauchy problem is via the 
method of characteristics. The characteristic lines x(t) are defined by the equation 



v t = a(v) v : 



v(x,0) = <p(x) , 



(6) 



x{t) 



a(v(x(t),t)) 



(7) 
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and any classical solution is constant along these lines, which are mutually distinct 
provided t* < t c , where 

t c = inf > . 

{a;| Vx a'( V )>0} Lp x a'(ip) 

The time t c is known as the critical time and at that time, the solution is said 
to have a gradient catastrophe: the solution remains bounded while its derivatives 
blow up. Being constant along the characteristic lines, the solution of the Cauchy 
problem ((6]) can locally be described by the functional equation 

x + ta(v)=g'(v), (8) 

where g is a function obtained by the implicit equation g'{tp(x)) = x. Equation l[8"]l 
admits a variational formulation, constructed in the following way: it is well known 
that for any function /, the (formal) functional H® = j f(u)dx is conserved along 
the flow of ©. In addition, one has that the family of functionals 



Q a , M := J (xa(v) +tf3(v) -/(«)) dx , 

is conserved, for any choice of /, provided the relation 

P'(v) = a'(v)a(v), 
is satisfied. Choosing / such that /' = a' g', then ([5]) is equivalent to 



8v 



0. 



provided a' does not vanish. Note that in this unperturbed case one can always 
choose a(v) = v, thus obtaining ([5]) directly. 

1.2 Hamiltonian perturbations 

Let h(u), c(u), p{u) and s(u) be smooth functions of one variable, with c not iden- 
tically zero. Then, introduce the quantities 

c h :=cti", 

s h := a ti" - — - C A h (6) _ £! ft (6) _ t h (A _ E h (5) 
8 8 24 6 6 

and consider the linear map defined by 

2 

D( c ^ s )(h) := h - e — c h (u) u 2 x + e 4 (p h {u) u 2 xx + s h (u) u*) , (9) 

where e is a real parameter. The map © has been introduced by Dubrovin in Q5], 
where he called it i D— operator'. Note that when applied to a quadratic function, 
^ reduces to the identity map. We therefore assume that h'" ^ 0, and we introduce 
the local functional 

H h := J D (c ^ s) (h)dx, (10) 
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the Hamiltonian of our system. The corresponding evolution equation, that is, the 
fifth order PDE given by 

' SHh 
Su 

is a Hamiltonian perturbation of the quasilinear PDE ^j. Explicitly, we have: 



u t = {u, H} = d x 



u t = h" u x + e 2 d x (J^c' h u 2 x + c h u xx ^j + e i d x (^2phU xxxx 

+ ip h u x u xxx + 2(p/ l - 6s h )u 2 x u xx + ip' h u 2 xx - 3s' h u{\ . (11) 

Proposition 1 (Dubrovin, [3]). The PDE pip admits infinitely many approximated 
conserved quantities (up to order e 6 ), of the form 

H f = J D M {f)dx, (12) 

parametrized by an arbitrary function of one variable. 

The above proposition can also be stated in the following way: for every /, g G 
C°°(K) the corresponding functionals Hf, H g satisfy: 

\ ou ax ou J 
where E is the Euler-Lagrange operator defined in ©. 
Example 1. The KdV equation, 

Uf u u x -\- s u XXX j 

can be obtained from the above class by choosing h(u) = c(u) = 1 and p{u) = 
s(u) = 0. Specializing to this case, we obtain that the generic 0(e ) — conserved 
quantity for KdV is given by (|12p with 

DkMI) = / - y f" ul + e 4 (± f^ ul x - i/W u% 

and f(u) an arbitrary function. 

The functions c, p and s appearing above characterize a hierarchy of Hamiltonian 
PDEs, which commute up to order e 6 . Every equation of the hierarchy is generated 
by an Hamiltonian (JTUJ) for a specific choice of h. As proved by Dubrovin in [1], 
the functions c and p (but not s) are invariants of the hierarchy with respect to 
Miura-type transformations of the form 

u^w = u + Y,£ k Fk(u;u x ,...,u^), (13) 

k>l 

with Fk polynomials in the derivatives of u, of total degree k. Equivalently, the class 
of Hamiltonian perturbations modulo Miura-type transformations (|13p is parametri- 
zed - up to order e A - by three arbitrary functions of one variable: h, c and p, while 
the function s can be chosen at convenience. In this paper, we present our the 
results for an arbitrary choice of the function s. Indeed, as shown for instance 
in the example below, many important examples of Hamiltonian perturbations are 
obtained selecting different choices of s. 
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Example 2. The following class of generalized KdV equations 

u t = a(u) u x + e 2 Ki u xxx + e 4 n 2 u xxxxx , 

where a is an arbitrary smooth function, k±, k 2 arbitrary constants, can be written 
as Hamiltonian perturbations by choosing 

, „ Kl K2 3 «? a" 

h = a, c= — , p — 



a' 2a' 10 (a') 3 ' 

2 (a") 3 7 a" a'" 1 a"" \ k 2 ((a") 2 



5 (a') 5 20 (a') 4 24 (a') 3 / 12 \ («') 3 («') 

1.3 Solutions to the perturbed equation for small times 

Let us consider now solutions of a perturbed equation ([lip . The discussion here will 
be at a formal level, for in general (that is, for an arbitrary perturbation) it is not 
known whether the equation admits a global - or even local - solution for a given 
initial data. We consider the Cauchy problem 

u t = {u,H h }, u\ t= o = <p, (14) 

where the Hamiltonian Hh is given by (|10p , and the initial data ip, which is e- 
independent, is the same as for the unperturbed equation A standard method 
to obtain approximate solutions to (|14p is by considering a semiclassical expansion 
in the small parameter e, which in our case reads as follows: consider the Ansatz 

u{x,t,e) = ^2v\x,t)e 2 \ (15) 

with coefficients v l smooth functions of x and t. Within this setting, all identities 
are understood in the sense of formal power series in e - they are assumed to hold 
identically at every order in e. In particular, by evaluating (|T5|) at t = 0, one gets 
the relations 

v°{x, 0) = <p(x), v l {x, 0) = 0, i > 1. 

By expanding both sides of equation (fTTj) according to the Ansatz (|15p . in first 
approximation one obtains 

v° t =h"(v°)v° x , u\ t=0 = V >, (16) 

which says that v° must be a solution of the Cauchy problem ©. Accordingly, from 
the higher order coefficients one obtains an infinite set of semilinear equations (or 
transport equations) for the coefficients v k {x,t). For instance, the equation for v 1 
turns out to be 

vl = d x (h"(v°)v 1 + c h (v°)v° xx +\c' h {v°) («2) 2 ) , t; 1 !^ = 0, (17) 
while the equation for v 2 is 
v 2 = d x (VV + ^'"(v 1 ) 2 + d x (c h vl) + c'.v^v 1 + \cl{vl) 2 v l + 2 Ph v° xxxx 

+ ¥ h v° x v° xxx + 2pUv° x ) 2 v xx + 3p' h (v° x ) 3 - 12s h (v° x ) 2 v° xx - 3^(«2)*) . (18) 
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with the initial condition v 2 \ t =o = 0. Here, the functions h, Ch, Ph, Sh and their 
derivatives are evaluated at u = v (x,t). Note that, although in principle every 
transport equation can be solved recursively starting from the solution of p^|) . the 
use of a standard method for finding the solution requires a case-by-case approach, 
which can be very difficult to implement, even in simple cases. 

Remark 1. It is reasonable to expect that expansion (|15p holds true in the space- 
time region where the solution of the Hopf equation Mb)) obtained by the method 
of characteristics is single-valued. This region can be described as the complement 
of the so-called 'Whitham zone' ]1 61 {!]/. For the class of perturbations of Example 
and for times smaller that the critical time, we recently proved in JWj that an 
expansion of the form (|15p actually exists at some order in e, provided the initial 
data is sufficiently regular. In the full complement of the Whitham zone, namely 
also for times t > t C) the existence of expansion (|15p was proved just for special 
classes of solutions of the KdV equation and its hierarchy - see J3\f for a review. 



An alternative way to study solutions of (|14p for small times has been proposed 
within the theory of Hamiltonian perturbations in [H [TT] . The idea is to look for 
a specific (quasi-)Miura transformation, mapping any solution of the unperturbed 
equation, say v°(x,t), to a solution u(x,t) of the perturbed equation. Remarkably, 
this is always possible for a e 4 — perturbation: the transformation is of the form 



v a .— ► u = v° 



+ e {v°{x),K}+e 2 {{v°(x),K},K} + ... (19) 



where the functional K, up to order 4 in e, is given by 

K = 



2 "* ~° " x \ 40 V v° J 4 v° 



dx. (20) 



Although this transformation turns out to be very useful when considering algebraic 
and geometric aspects of the theory of Hamiltonian perturbations [H [10], if one 
is interested in (approximate) solutions of (jlip . then the above formula has some 
drawback. Indeed, it is clear that (fT9"ll2T))) is singular at critical points of solutions to 
the equation ([5]); moreover, the required initial condition Uu—q = tp is not satisfied. 

In order to eliminate these problems, we construct in the next section a suitable 
deformation, to the case of Hamiltonian perturbations, of the characteristic equation 
l[5]). and we show how to use it in order to explicitly compute the correct terms of 
the semiclassical expansion. 



2 The variational string equation 

The aim of the present section is to construct an analogue of equation |S} for Hamil- 
tonian perturbations of type (fTTj) . The main idea we follow is to suitably generalize 
the variational formulation of the method of characteristics, outlined in Section [TJ to 
the perturbed case. Remarkably, it turns out that this is always possible at order e 2 , 
while at order e 4 nontrivial constraints appear. There are three properties that we 
require in order to obtain a meaningful deformation: first, the new equation has to 
be written as the critical point of an (almost) conserved quantity for the perturbed 
equation. Then, at e — we want to recover the method of characteristics, and 
finally, in the limit as t tends to zero we want to get the correct initial data. If (p 



S 



is the initial data of (jlip . g a function such that g'(ip(x)) = x, and S denotes the 
functional we want to find, then these conditions can be written as 

d t (^) + £ ( IT T"^) = °{z K l fOT some K > 0, (21a) 
\ou J \ ou ax ou J 

^ =0 x + h"(u)t = g'{u), (21b) 

OU | e =0 

^ =0. (21c) 

OU \t=Q,u=ip(x) 

Maybe only equation (|21ap requires some comment. If S is an exact conserved 
quantity, explicitly depending on t, for equation pip , then the right hand side 
of (|21ap is exactly equal to zero. Since in this case we expect to obtain only an 
approximate conserved quantity, then (|21a[) tells us that the critical locus of S 
varies, at least formally, slowly: we can expect the solution of the Cauchy problem 
to satisfy the critical point equation up to an error 0(e K ), uniformly in x, t on some 
space-time domain. 

Given a Cauchy problem for a perturbation (jlip , if a special conserved quantity 
S satisfying conditions (|2ip is found, at some order in e, then the required analogue 
of the characteristic equation © is given by 

f=0. (22) 
ou 

We call this equation variational string equation. In the special case when the 
perturbed equation (jlip is given by the choice c(u) = co,p(u) = po, s(u) = 0, with 
Co and po constants, a possible candidate appears in the Dubrovin's paper [3], where 
he considered the functional 

S = J (xu + tD Co , Pofi {ti) - D COiPO! o(gfjdx, (23) 

in relation with the universality conjecture of solutions at the critical point. Note, 
however, that the above functional satisfies condition (|21cj) at order 0(e 2 ) only. 

In the following, we generalize the functional (|2"3"|) to a more general class of per- 
turbations. Note that this quantity is the sum of two separate conserved quantities: 

S = J (xu + tD comfi (h'f)dx, (24) 

and 



H g = J D COtPofi (g)dx. 



We will first consider an extension of the former, requiring only (|21ap and (|21b[) 
to be fulfilled. An appropriate choice of the second conserved quantity will then 
provide the correct initial data (|21cp . 

Remark 2. The reason of the name 'variational string equation' is due to the fact 
that in a particular case identity (|22p reduces to the 'string equation' associated to 
the KdV hierarchy ] 151 \13jf . Let us consider the critical point of the following KdV 
exact conserved quantity 

S = I [ xu + t— J dx — H n , 
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where H n is the n-th Hamiltonian of KdV hierarchy (the deformation of J u n dx, 
conserved quantity of the Hopf equation). The variational string equation for this 
functional reads 

, . SH n 
x + ut = — — . 
ou 

Note that this is an exact identity. Differentiating the last equation with respect to 
x, we get 

l = -tu x + ± ^ = -[L,P], P = -{L^) + + t 
ax ou 

where L = —£ 2 -§^i +u is the Lax operator of the KdV hierarchy, and {L^~) + is the 
differential part of the pseudo-differential operator L~^~ . The equation [P, L] = 1 
(after a suitable rescaling x — ?> ^) is known in the literature as string equation. 

2.1 A special conserved quantity 

Consider a family of functionals, depending on t, of the form 

(xp(u;u x ,u xx ,...,e) + t p,(u;u x ,u xx , . . . ,e)J dx, (25) 
for certain differential polynomials 

p =p + e 2 (piu 2 x + p 2 u xx ) + e 4 (pzu x + piU 2 xx + p 5 u 2 x u xx + p e u x u xxx + p 7 u xxxx ) , 
p =p a + e 2 (piu 2 x + p 2 u xx ) + £ 4 (p3U x + PiU 2 x u xx + p 5 u 2 xx + p 6 u x u xxx + p 7 u xxxx ) , 

where the p, and pi axe arbitrary functions of u. We want to determine conditions 
on these functions in order to obtain a conserved quantity - up to order 0(e 6 ) - for 
equation pip . As a first step, we write this functional in normal form, that is up to 
exact derivatives: 

Lemma 1. Every functional of type (|25p can be written in a unique way in the form 
J (x a(u; u x , . . . , e) + e 4 e(u) u x + t f3(u\ u x , . . . , e)^ dx, (26) 

where 

a =a + £ aiu x + e (a 2 u x + a 3 u xx j , 
p =Po + £ 2 Pi u 2 x + e 4 (fa ui + f3 3 u 2 xx ) , 

and oti(u), (3i(u) and e{u) are functions of one variable. 

Proof. A direct computation, using integration by parts. □ 

Proposition 2. For arbitrary functions h,c,p and s, with c not identically zero, 
equation (jlip admits a conserved quantity of order 0{e 2 ) of type (|26p . This is 
specified by choosing a = D^ c p ^{r) and (3 = D^ c ^ p s j(q) in (|26p . where the functions 
r(u) and q(u) satisfy 

1 h" 
r' = ~, q' = ~. (27) 
c c 

The 0(e 2 )- conserved quantity is unique up to a multiplicative constant. 
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Proof. A direct calculation: consider a functional of type ([2"5[) and impose that it 
is conserved at order 0(e 2 ) with respect to equation (|TTj) . Denote uq = r, 0q = q. 
Then, at order O(l), one gets the relation 

r' = q'h", (28) 

while at order 0(e 2 ) there are three conditions: the first two, 

1 III a III 

u\ = --cr , Pi = --cq , 

imply that - up to order e 2 - we have a = Dr CtP S \{r) and f3 = -D( c ,p,s)(?)- The last 
condition is given by 

3r"ch"' + r"'ch" + r'ch}- 4 ' + r'ch'" — cq'" = 0, 

and after substituting ([28)1 into it, one obtains (|2"7)l . □ 

Let us consider now the 0(e 4 ) term. 

Proposition 3. For arbitrary functions h, c and s, with c not identically zero, 
equation ([lip admits a conserved quantity of order 0(e 4 ) of type ([26p if and only if 

p(u) = lc(u)c'(u) + Xc(u) 3 , (29) 

5 

where X is an arbitrary constant. In this case, we have 

_53 c"{c'f 8 (c'f £^157 (O! _ , _ 173 cV^ _5_ (4) 
S ~24 c 2 15 c 3 + c 120 c c 2 120 c + 12° 

, / 13 „ . 2 (c') 3 11 „A 
+ A^ y cV + ^ + - cc '"j. (30) 

T/ie 0(e 4 )- conserved quantity of type ()26p is given as in Proposition^ and it is 
unique up to a multiplicative constant. 

Proof. As for Proposition^ the condition of ([2U)) being a conserved quantity at order 
£ 4 implies that a and /? are themselves 0(e 4 )-conserved quantities. In addition, one 
finds two more conditions. One of them involves the function e(u) and is satisfied 
by the choice (|30p. the other reads 



h^_ 



(p'c-'3pc' + ^c(c') 2 -^c 2 c'^ =0, 



from which one obtains relation ([2"9"[) . □ 
Due to the above propositions, the required functional has the form 

So = J (xD C:P , s {r)+tD (c ^ s) (q)+e 4 e{u)ul)dx, (31) 

with r and q given by ([2T[) and e given by ([3"0[) . By construction, condition (j21a[) is 
satisfied for any Hamiltonian perturbation at order e 2 , and for the class of equations 
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characterized by (j2"9")l at order e 4 . Note that since the density of So depends explicitly 
on x, its variational derivative is given by 

^ = X £(D c , pA r)) + + E(tD (CtP , s) ( q) + e*e(u)ul). 

k>0 

Example 3. In the particular case c(u) = co,p(u) = po,s(u) = 0, then the func- 
tional (|3ip reduces to (|24p . Note that (|29p is satisfied by choosing A = po/cg- 

Remark 3. One of the hypothesis on the Hamiltonian perturbation in Proposition^ 
and\5\is that c(u) does not vanish identically. Indeed, if c = then the proof breaks 
since the special conserved quantity So is built on the function -. Although we do 
not treat the case c = in generality, the following simple example shows that the 
string equation can be considered also in this case. The generalized KdV equation 

obtained choosing c = 0,p = |, s = in (|10p . admits the exact conserved quantity 

S = J |m + fyj dx . (32) 

Note that Sq is nothing but the conserved quantity associated to the Galilean invari- 
ance of the equation. 

Remark 4. A natural question arising in the above construction is whether the 
constraint (|29p at order e can be eliminated by a different choice of the class of 
conserved quantities (e.g., depending rationally on the derivatives). This, however, 
is outside the scope of the present paper. 



2.2 Fixing the initial data 

We now address the problem of recovering the required initial data, namely to 
modify (|3"Tj) to get a functional which satisfies (|21cp . Consider a Cauchy problem 
for equation (fTTjl . with given initial data u(x,0,e) = <p(x), independent of e. In 
addition, we introduce the function /, defined by the relation g' = c f . Due to 
the definition of g, we thus have that / is locally related to the initial data by the 
identity 

c(<p(x))f'(<p(x))=x, (33) 

for x belonging to some interval of the real axis. The correction of the initial data 
(in this case, at order at most e 4 ) can be obtained as follows: for arbitrary functions 
/i and /2, we consider the functional 

H f+£ 2 fl+E 4 h = J I ),,,,. J f + e 2 fi + e 4 f 2 )dx, 

and then truncate it at order e 4 , by introducing the 4 order (in e) polynomial Hj, 
defined by the relation Hj — Hf +e 2f 1+E if 2 = 0(e 6 ). Finally, define 

Sf.= So-Hf, (34) 
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which by construction satisfies (|21a[ 1 and, due to the definition of /, also (|21b[) . Note 
that (|3"4")l itself is a polynomial of order 4 in e. A suitable choice of the functions f± 
and fa implies condition ()21c[) . 

Proposition 4. For any Hamiltonian perturbation of type ill}) with arbitrary coef- 
ficients h, c 7 p, s the functional (|34p with So given by 1125)). is conserved up to 0(e 4 ). 
Moreover, if f is chosen such that (|33|) holds and 



c(v) 



2c(v) dv 2 \(c{v)f'(v)y J ' 



(35) 



then we have 



5u \t=0,u=<p(x) 



0(8*). 



Proof. The fact that the functional (f3"4")) is a 0(e 2 )-conserved quantity of (fTTj) for any 
choice of the functions / and f\ follows immediately from the definition of the linear 
operator _D c PjS . Let us prove the second part of the proposition. The condition 



5u \t=0,u=tp(x) 



0(e 



is verified due to the choice (|33p . Moreover, the e 2 coefficient of Sf is given by 



/ 



h + o [x - 



dx. 



The variational derivative of this functional, evaluated at t — 0, u — tp(x) and set 
equal to zero gives the following condition on f\\ 



2(d 



>V2 



C 

+ - 7) Vx(s) - c/"Vxx(x) + i (c'/'" + c/' 4 ') ^(*) 2 . (36) 

Of course here the functions /, c and their derivatives are computed at u = (p(x). 
Using condition (|3"3")) . together with the differential identities 



5 c'c" 



2 c 



1 



(c/0 



,y |u=<p(x)> 



%x(a;) 



obtained from it, one gets that ([31))) is satisfied by the choice (|3"5|) . In particular, 
this implies 

OM \t=0,u=<p(x) 



and the thesis is proved. 



□ 



2 at order e 2 . The identity at e 4 is verified subject to condition (1291) . 
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Proposition 5. Under the hypotheses of the previous proposition, and if in addition 
the coefficient p of the Hamiltonian satisfies condition (|29p . then the functional (|34l) 
is conserved up to 0(e 6 ). Moreover, if fi is chosen such that 



i d 2 

c{v) dv 2 



h(v) 



whe 



f - P 



[(cf')fdv 2 \(cf>)' ) dv\[( c f>) 



J d 



[(c/01 



1 



( / :s 



c c 



6[( c /')f dv* \{cf>yj 24[( c //)f cfo 2 V(c/')' 



i^ cc ' c "/'J_ c 3 / (4) 

120 J 10 J 



J_ c V"/'--(c') 3 /' + — c 2 c7'" 
10 5 V ; ^ 60 ^ 



30 v 7 J 8 J 



1 



[( c /')fdw V(c/')' 



1 



+ 



1 / /\2 /■/// i 1 / // pi I 1 // till 

o(c) / +2 C c / ~4 CC / -~ cc f + 7^ cc / 



1 / 
12 C 



i c 2 /(S).!^) J/ 

24 ^ 24 ^ 



[(c/')'] 



<5S 



= 0. 



then we have 

5u \t=o,u=ip(x) 

Proof. The proof is similar to the previous proposition. 



□ 



The functional (|34[) with f\ and f2 given as in the above propositions is the 
required functional for the variational string equation, as it satisfies conditions (|21|) 
at the desired order. The variational string equation (of order K) associated to Sf 
is defined as the equation for the critical point: 



5S 



f 



0. 



(37) 



Explicitly, at order 0(e 2 ), the variational string equation reads 
1 



(x + th" -cf)+s 2 { x 



>=0(e i ). 



+ t 



h"\" 



h"\"\> 



It is a nontrivial and remarkable fact that this equation provides approximate solu- 
tions to a nonlinear PDE for arbitrary initial data. 
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It is natural to ask in which case the construction of the variational string equa- 
tion can be extended higher orders, and more generally how to select Hamiltonian 
perturbations admitting a variational string equation at any order. There are two 
main difficulties to solve this problem: first of all, for a generic perturbation and 
for an arbitrary function /, it is not known whether the functional J f(u)dx can be 
deformed at any order to a conserved quantity with density depending polynomially 
on the jet variables, as required in the construction of the variational string equa- 
tion. In the language of |11) . this requirement is the definition of formal integrability 
for the perturbed equation. In addition, it is in general not known whether special 
conserved quantities of type 5*0 exist at any order. A positive answer to these prob- 
lems is given by the KdV example. Indeed, due to Galilean invariance, KdV admits 
the exact conserved quantity ([52"]) , and the equation itself has been proved in [TT] to 
be formally integrable. Therefore, the string equation in this case can be extended 
at any order. 



3 Approximate solutions to perturbed equations 

We now use the variational string equation (|37p to compute the first two perturbative 
corrections of the solution u(x,t,e) to an Hamiltonian perturbation (fTTjl . In |12) . 
we gave the first correction v 1 without explaining in detail the procedure followed 
to obtain it. Here we fill this gap, and compute explicitly the second term v 2 . 

Suppose the Cauchy problem of (fTTj) admits the semiclassical expansion (fT5|) . 
and suppose in addition that the equation admits a variational string equation up 
to order 0{e 2N ). Then, we first require u(x, t, e) to be an (approximate) solution of 
the variational string equation: 



Su \u=u(x,t,e) 



0(e 2N+2 ), 



and we expand the above formal identity according to the Ansatz (|15p . At order 
0(1), we find 

x + th"(v°) = c(v°)f'(v°), (38) 

which is precisely the formula for the method of characteristics, thus implying that 
v is a solution of the dispersionless equation, as required. Consider now the e 2 
term. This reads 



xc' + t(h"c' -h'"c)+c 2 f" 



■ I - 



,,() 



: I - 

c 



(v° x ) 2 



Q\2 



■■f"v xx + (cf")'^Y 



1 



(cf)' 



+ t 



h"\» 



h'\»y(v° x y 



where all functions are evaluated at u — v°(x,t). The explicit dependence on / of 
the above formula can easily be eliminated: further differentiating the exact identity 
(|38p with respect to x, one obtains a formula for the quantities f^'(v°). Substituting 
into the equation for v 1 , this gives 



v 1 (x,t) 



t d f {ch'")' (v° x ) 2 + 2ch 
2 dx\ 



/// o 



+ tc(h'") 2 v x v xx +tc'(h"r(v x ) 3 



(l + th"'v° x ) 2 



(39) 
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which is a universal formula: given the solution to the dispersionless Cauchy problem 
(|16p. then (|3"9"f satisfies the transport equation (|17|) with the correct initial data. 
Note that using the variational string equation the problem is solved by a finite 
number of elementary manipulations. A more elegant representation of v , first 
introduced in |12| , is given by considering the functional (depending on t) : 

Kx[u] := -i J c{u) u x log (l + tti" (u) u x ^ dx. (40) 

Then, the function v admits the Hamiltonian representation 

..(,,», »(«jm y 

dx y 8u{x) lu=v0{x ^j 

The calculation of the higher order corrections can be computed similarly, provided 
the variational string equation associated with the PDE exists at the desired order. 
For instance, let us consider the second correction term for the Hamiltonian satisfy- 
ing condition (|29[) . Since the formulae in this case are much more cumbersome, we 
present here only the more compact variational formulation. Indeed, it turns out 
that the transformation from v°(x, t) — > u(x,t,e) can be written - up to order e 4 - 
in the canonical form 

u = v° + e {v°, K) + ^{{v°, K], K } + . . . , (41) 
K = sK x + s 3 K 3 + . . . , 
where Ki is given by (|40p and K% is the following expression: 

KM = _( ( c V" + c "' \ cu x log(l + th"'u x ) 



+ 



{h'") 2 ) {h'") 2 J 8t 2 

1 c 2 {ti")H 3 ul x (5 + 5 1 ti" u x + t 2 {h!"f u 2 x ) 
40 (1 + th"'u x ) 5 

3 c 2 t h^u 2 xx (4 + 10 1 ti"u x + 5 1 2 (h"') 2 u 2 x ) 1 th'"su x {2 + th"'u x ) 
40 (1 + th"'u x f + 2 (1 +th"'u x ) 2 

1 (3cc' + 5Ac 3 ) t h'"ul x {2 + tti"u x ) (2 + 2tti"u x + t 2 {h"') 2 u 2 x ) 



4 



20 (l + th"'u x ) 

)^ c2 \ (2ch^ +&dh^ 3tc{h^) 2 u 



240 1 



20 12 J V {l + th"'u x f (1 + th'"u x ) 4 i 

12 c 2 (/iW) 3 3ch^ (3 h'" h^c' + 3 h^h"'c- 16 c(M 4 )) 2 ) 

{h'") 4 {l+th'"u x f + (h'") 4 {l + th"'u x ) 4 

3hW (8cc"(h"') 2 -4c 2 (h^) 2 + 9cc'h'"h^ +9c 2 h'" h& - 16 (c') 2 (h'") 2 ) 

(h'") 4 (i + th"i u x y 

96 {c') 2 {h"'fh^ - A8cc"{h"') 2 h^ + 63cc'/i'"(/i (4) ) 2 + 5cc"'(ti") 3 

(h"') 4 (l +th"' u x ) 2 
252c 2 (feW) 3 -GZWhWh* -30 c 2 (h"') 2 h^ + 441 c 2 h'" 
{h m Y{l + th'" u x ) 2 + (h"') A {l +th"' u x ) 
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+ 



90 cc'(h"') 2 h ( V - Mlcc'ti" (h^) 2 + 768 c 2 {h^f + lUcc" (h"') 2 h^ 

(h"')\l+th"' 11 x) 

48 {c') 2 {ti") 2 h^ +35cc"'(h"') 3 ISOcfeW / (4 ,y 90c / , U) 
(h"') 4 (l + th"'u x ) + {h"'Y V J {h'") 2 

30 c 2 ^ 6 ' 180c 2 (/iW) 3 35 cc'" 5cc' 



(h'") 2 



-(1 + tti" u x ) 



The canonical transformation (|4H . with K\ and K$ as above, is well-defined for the 
solution v° of the unperturbed equation up to the time of the gradient catastrophe. 
The apparent singularities of K3 at t = and h'"(u) ~ cancel out; in fact, K3, 
as well as K±, vanishes identically if t = or h'" = 0. In the particular case when 
h"{u) = u, the above formula reduces to 



K 3 [u] 



1 cc"'u x log (1 + tu x ) 1 t 3 c 2 (u xx ) 3 (5 + 5tu x + t 2 (u x ) 2 ) 
4 1? 40 (1 + tu x f 

1 tc(3c' + 5c 2 X)(u xx ) 2 {2 + tu x )(2 + 2tu x +t 2 (u x ) 2 ) 
+ 20 (1 +tu x ) 4 

1 (2 + tu x ) [cc"'(6 + 6tu x -t 2 {u x ) 2 ) +24st(u x ) 2 ] 
+ 48 t(l + tu x ) 2 ' 

A direct computation shows that the corrections v , v 2 computed from the canonical 
transformation (|4ip satisfy the appropriate transport equations (|171 1 . This is an 
a posteriori confirmation of the validity of the variational string equation. 

Example 4. For t 6 [0,t c ), the solution to KdV with initial data (p is given, at 
order e 4 , by 



u = v u + -dUlog(l + tv») 



2 

td. 



(l + t«2) 
3 20 + 15tv a x + 3t 2 (v°) 2 



e 4 t 2 



t 2 (5 + t^)(^ ) 3 
10 {l + tvlf 



10 



{1 + tvlf 



0. 



+ n 

x \{i + tvi) 4 x 



+ 0(e e ), 



where v° is the solution of Hopf with same initial data. Note that att = 0, the above 
formula gives u(x,0,e) = v°(x,0) = <p(x). 



4 The string equation is an approximate identity 

In this Section we discuss in which precise sense solutions of a Hamiltonian pertur- 
bation of the Hopf equation satisfy a variational string equation. Indeed, given the 
Cauchy problem (fT4")l and the functional Sf conserved up to 0(e K ), we will prove 
that the solution u(x, t, e) of the Cauchy problem approximatively satisfies the string 
equation: 

6 -^(u(x,t,e)) = 0(s K ), 

provided it is sufficiently regular in the e = limit. The problem presents many 
difficulties, as the string equation is the critical point equation for an ill-defined 
functional (in the generic case, g'(0) 7^ 0, the integral diverges for any rapidly 
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decreasing function), which is not exactly conserved. Moreover, the string equation 
itself is singular and valid only locally in the space-time plane. 

The first step of our proof is to show that the variational derivatives of a formal 
functional satisfies a linear PDE with a forcing term, which vanishes in case the 
functional is (formally) conserved. In the case of the string equation, this forcing 
term will be small, proportional to e K for some K > 1. 

After this result, we are lead to analyze local solutions to linear PDEs with 
a small forcing term. Here, the standard methods of functional analysis are not 
effective, as the equation is valid on some bounded domain of the real line and we 
ignore the boundary values. Equivalently, when estimating the time derivative of 
the L 2 norm of the solution and integrating by parts some unknown boundary terms 
appear. 

To overcome this problem, we will estimate the (time derivative of the) L 2 norm 
on an interval whose extremes vary with time. Remarkably, if the extremes vary 
with a simple law - reminiscent of the characteristics' equation - the boundary terms 
cancel-out. 



4.1 A linear equation for conserved quantities 

Here we prove that the variational derivative of a formal functional M satisfies a 
linear equation with a forcing term which vanishes in case the functional is formally 
conserved. 

Theorem 1. Let M be a formal local functional, possibly depending explicitly on 
x,t, let fh = be its variational derivative, and let u(x,t) be a solution of the 
Hamiltonian PDE 

d SH 

dx Su 

Then, the function of two variables Ai{x, t) := fh(x, t, u{x, t), u x (x, t), . . . ), satisfies 
the linear equation 

— = J2A l (x,t)^- Tr +K(x,t), (42) 

i>0 



whe 



re 



. , , ( d SH\ „, . ( n SM ^,SM d SH, 

provided all the partial derivatives on the right hand side exist and are continuous. 
Proof. For any sufficiently regular solution u, we have that 



~dT~\ " + dW)~6u ~dT 



\u— u(x,t) 



SM d SM d l+1 SH 

1 Su Z— 1 9u(») Su dx l+1 Su 

i>0 
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The thesis follows from the remarkable identity due to Gel'fand and Dikii [TJ: 

d f g+l gg - -rf 

^du&dx i + l9 2-<duWdx i + lJ KJ dx 9h 

i>0 i>Q 

which is valid for any two /, g variational derivatives of formal local fvmctionals 

j — su ■ y — su- lj 

Let us analyze the result of the previous theorem. In case M is formally con- 
served, equation holds with 1Z = and it therefore reduces to a linear PDE. 
In this particular case, the equation was already obtained by Lax in [9]. Theorem 
Q] above extends Lax' result in two directions: first, we do not assume M to be con- 
served, secondly, the equation remains valid in case M is just a formal functional. 
Such a generalization is important for us, as the variational string equation is the 
critical point equation of a formal functional Sf that is not strictly conserved. 



4.2 The domain of the variational string equation 

We begin now to give a precise sense to the variational string equation; as a first step, 
we define the domain where the equation is defined. Let us choose an Hamiltonian 
perturbation of type (jllj) . admitting a string equation up to a certain order in e. In 
addition, we fix: 

• A solution u{x, t; e) of the corresponding Cauchy problem 

u t = d x -r± = h"{u)u x + 0{e 2 ), u(x,t = 0,e) = ip(x). (43) 
Su 

• An interval (x—,x+) of strict monotonicity of ip such that c(cp(x)) ^ for all 

x e {x-,x + ). 

Let /' : (a, b) — > K, be defined by the relation 

f'(<p(x)) = , x£(x-,x+). 

and let 5/ be the special quasi-conserved functional constructed in Section [21 We 
recall that Sf satisfies the properties ([21"]). which can now be written as: 

dt —r^- + £ ( —r-^- — — - — )=e + r , for some N > 1. (44a) 
ou ou dx ou 

5 -§l =J-( x + h "(u)t-c(u)f(u)) (44b) 

OU |e=0 C(U) 

SSf 

=0, x£(x-,x + ) (44c) 

OU \t=0,u=ip(x) 

Here r is some differential expression, polynomial in e 2 . We can now start to describe 
the domain of the string equation. Let us denote by a, the variational derivative of 
Sf evaluated at u = u(x, t; e), that is: 

SSf 

a(x,t,e) := —^-{u{x,t;e)) . 
Su 

Then, the following facts hold: 
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i) The function o~(x, t, e) is well defined whenever u(x, t, e) belongs to the domain 
of /', namely the interval (a, b). Indeed, the only source of singularities is the 
zero locus of c(u) and, by construction, c(u) ^ on (a, b). 

ii) By construction, a(x,t,e = 0) = on the domain of the [x, t) plane filled by 
the characteristics line emanating from the interval £+). Indeed, on that 
domain equation ([SJ of the method of characteristics holds. 

iii) Due to (]44c(l . we have a(x,t = 0,e) = if x £ 

A meaningful domain for a must contain the subsets described in ii) and iii) , where 
er is known to vanish. In addition, we ask the domain of a to be connected. 

Definition 1. Let f2 = {(x, t, e) € K 3 | \u(x, t, e) £ (a, b)} and let f2 c be the connected 
component of fi containing (x—,x+) x {0} x R. We define f2 c to be the domain of 
a, and we denote by f2° o the intersection of f2 c with the hyperplane e = Eq. 

Notice that fi§ is exactly the domain filled by the characteristic lines (|7|) of the 
unperturbed equation with initial point xq £ (x_,x+). 



4.3 Proof of the validity of the variational string equation 

In Theorem [I] we have shown that the variational derivative of a formal functional 
satisfies the linear PDE (|42p with a forcing term proportional to the total time- 
derivative of the functional. Combining this result with (|44a[) . we obtain that on 
the domain fi c , the function a satisfies 

E^S +£2W+2 ^ *<=° = ' ( 45 ) 



dt ^ dx l+1 

i>0 



where 



\6uW 5u J [u=u(Xtt . £) 

and r is as in (|44ap . Notice that r is a polynomial in e, so that the forcing term in 
(|45p is 0(e 2N+2 ) provided u is continuous as e — > 0. If we consider just the explicit 

dependence on e 2 , then the linear differential operator <Ai g x t+i acting on a is a 
polynomial in e 2 . In particular, we can write it as: 



d l+1 ,„, x d 



V\ /1 

i=0 

where h"(u) is the same as in and A is a linear differential operator polynomial 
in e 2 , depending on u and on its x-derivatives. Using (|45[) and above decomposition 
(1461) . we prove that the string equation is satisfied up to 0(e 2N+2 ). 

Theorem 2. Suppose J^5[ ) depends on u as well as on its first L derivatives, and 
let [0, T] be a time interval such that the quantities 

di+ 2k u(x, t, e) 
dx^de 2k 

exists and are continuous for j < L, k < N for (x,t,e) £ f2 c . Then, we have that 
a = 0(e 2N+2 ) on compact subsets of £l c n {t £ [0,T]}. Therefore, the conditions 

d 2k cr(x,t,e) =Q 



ds 



2k 



=0 
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hold for k = 0, . . . , N and for any t G [0, T] . 

Before proving the theorem we need a lemma. 

Lemma 2. Under the same hypothesis as in Theorem^ let x(t;e) be the solution 
of the Cauchy ■problem 

dx 

— = -h"(u(x(t),t;e)), x(t = 0, e) = x(0) G x+), (47) 

with h" (u) as in |^3| ) and initial data x(0) independent of e. Then, we have: 

(i) For any initial data x(0), the point (x(t; e), t, s) belongs to Q c £ for any t G [0, T], 
provided e is small enough. 

(ii) Letx\,xi be two solutions of J^?] ) with initial data Xi(0) < £2(0), thenx\(t\e) < 
X2(t;e),Vt G [0, T] provided e is small enough. 

(Hi) Let K be a compact subset of £l c and K £ * := K n {t G [0,T],e 2 < e' 2 }. Then, 
there exist Xi(0), £2(0), with X\(0) < 2:2(0), and e' > such that 

K e ,C (J ([x 1 (t,e),x 2 (t,e)],t,e) C f! c . 

e 2 <e' 2 ,t£[0 : T] 

.ffere [2:1 (i,e), 22(4, e)] is the closed interval with extremes xi(t,£),X2(t,e). 

Proof. We prove (i). Points (ii,iii) can be proven along the same guidelines. Along 
the curve x = —h"(u(x(t),t;e)), the function u varies slowly, namely 

— = e 2 P(u, u x , . . . , , e) , for some j < L , 

where P is a differential expression , smooth in u and its first j derivatives and 
polynomial in e 2 . Therefore u(x(t, s), t, e) — u(x(0),0,e) = 0(e 2 t) uniformly in 
t G [0, T]. Hence if e is small enough, u(x(t, s), t, e) belongs to (a,b), the domain of 
/'. This proves the thesis. □ 

Proof of Theorem [H The principle of the proof is to show that the L 2 -norm of a is 
0(e 2N+2 ) (uniformly) when evaluated on any subinterval of ([xi(t, e), x^{t, e)]), for 
some function x\(t,e), and X2(t,e). Let us compute how the L 2 -norm varies on a 
time-dependent interval. Using (|45l l46|) . we get: 



dt 



/xi(t) a2dX \ = /*l 2(7 ( h "( U )°* + e2ACT + e 2N+2 TZ)dx + X 2 a 2 (x 2 ) - X\CT 2 (^1 ) 



= I _ Q h"'(u)u x a 2 + 2a(e 2 Aa + e 2N+2 K))dx 

+ [x 2 + h"(u(x 2 ))] a 2 (x 2 ) - [xi + h"(u(xi))} v 2 (x!)\ 
< sup xe[xi X2 ] \h"(u)u x \ f^ cr 2 dx + \x2 + h"(u(x 2 ))\ cr 2 (x 2 ) + \x\ + h"(u(xi))\ cr 2 (xi)+ 



Here Xi = "X ' and the operator A is defined in (|4"fJ)) . Choosing the boundary of 
the interval to evolve according to the law (|4"T)l . we get rid of the a 2 ( X i) terms which 
are not controlled by the L 2 norm[§ 



ft a2dx ^ su Px€t*i,« a ] ft "( u ) Kl C (72da; + 

(£ . 2 ^) 1/2 ( £ 2 (C(A CT ) 2 dx) 1/2 + £ 2W+2 (C fc a «fa) 



1/2 



3 Note that in the e = 0, this is the law of the characteristic lines Q. 
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We define 



n 



(J ([x 1 (t,e),x 2 (t,e)],t,e) 

e 2 <S 2 : te[0,T] 



for some 5 > 0. Applying the Gronwall inequality (recalling that <ti 4 _o = identi- 
cally), we get the following bound for a: 



x 2 (t) 



a dx 



Xl(t) 



< e 2 A(B 



,2N 



DT 



D 



(48) 



A = sup \x — y \ 2 < oo , B = sup |Ao"(x, t)\ < oo , 
x,yeu n 

C = sup \1Z\ < oo , D = sup \ h"(u)u x \ < oo . 
n n 

Since by Lemma [2] (ii) we have that x\(t) < X2(t) for any t, provided e is small 
enough, the above estimate implies that a = 0(e 2 ) on II. Therefore, we have 
a = 0(e 2 ) on any compact subset K of il c . Indeed, the previous estimate holds for 
arbitrary xi(0), £2(0) and we know, by Lemma [5] (iii), that K n {e 2 < e' 2 } C II for 
some xi(0), 2:2(0), and e' . Therefore, by continuity, a = on fig. We can now prove 
the thesis by induction. Suppose 



a 2 v 



0. 



j = 0,...,k 



de 2 i |e=o 

for some k < N — 1. By the hypothesis on differentiability, 



0. 



e=Q 



J = 0, . 



and so Act = 0(e 2k+2 ) on compacts of fi c . Hence, ((48]) holds with B = B'e 2k+2 for 
some constant _B'. Therefore, 



x 2 (t) 



xi(t) 



r 2 dx 



0(e 2k+i ) 



for any t £ [0,T]. Reasoning as above, we obtain 

d 2k+2 a 



Q £ 2k+2 



0. 



The induction step is proven. 



□ 



Remark 5. In Theorem&we can relax hypothesis h44 c \ l- Indeed the same Theorem 
holds if we choose an initial data if such that ^f { _ u _ = 0(e 2N+2 ) . The proof 
is valid with the only immaterial modification of adding an 0{e 2N+2 ) term on the 
right hand side of the estimate 
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5 Particular examples and quasi-triviality for KdV 



We have just shown that the String equation is approximately satisfied, i.e. up 
to 0(e 2N+2 ) - provided the solution of the Cauchy problem is regular enough in 
the semiclassical limit. In fact, the hypothesis of Theorem (J3J) is rather reasonable. 
We expect it to hold for very general Hamiltonian perturbations of Hopf with a 
sufficiently regular initial data and for any time smaller than the critical time t c . 
However, the only rigorous results in this direction are in our previous paper [12] . 
where we consider the Generalised KdV case, 

M 

u t = h"(u)u x + £ 2J a 3 u (2j+1 \ Oj G R, (49) 

u(t = 0) = ip e H S (R), s>2M + l. 

Here H s is the Sobolev space of index set. We proved that the Cauchy problem 
has a well-posed semiclassical limit, namely 

Theorem 3. Fix K G N, K > 1 and let u(x,t,e) be the solution of the Cauchy 
problem If the initial data ip belongs to H S (M.) with s > 3MK + 2M + I, 

then, for a time T > small enough, the following partial derivatives exist and are 
continuous for every x,t,s, with t G [0, T]: 

d i+j u(x,t,e) 1 

„ V .' ' , i<s-3MK- -, j<2MK. 

In particular, if ip G ris>o-^ S (^) an d f or ( x ^) S R X [0, T), then the solution u 
admits an asymptotic expansion as e — » 0, where all the coefficients are smooth 
functions of x andt. 

This is a corollary of the more general Theorem 9 in [T^]. We can apply this 
Theorem to prove that the string equations is indeed satisfied - up to the given 
error - by the solution the generalized KdV (|4U)) . From the theory developed in the 
previous section we can distinguish three cases: 

A Any generalised KdV equation admits a string equation of order 0(e 2 ). 

B In some cases, for example if h is cubic in u, then the generalised KdV equations 
admits - at least - a variational string equation of order e 4 . 

C Some particular equations, namely KdV or mKdV, which admit a variational 
string equation of arbitrary order. 

In all cases, combining Theorems [5] and [3] we can draw the conclusion that the string 
equation is indeed satisfied up to the corresponding error, provided the initial data 
of the equation belongs to a Sobolev space with high enough index. To give an 
example of how high the index must be, we consider the string equation of order £ 4 
for KdV and mKdV. In this case, the linear equation (|4"5) l satisfied by a depends on 
u and its derivatives up to Due to Theorem^ the String equation is satisfied 
up to 0(e e ) if 

d 6 d 4 u 
dx 6 ds 4 

exists and is continuous. From Theorem [21 this holds provided s > 
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Remark 6. Suppose the Cauchy problem J^ff| ) is globally well-posed if e 7^ and 
e small enough. Then the time T of Theorem [5| can be chosen to be any time 
smaller than the critical time of the solution of the Hopf equation. The above situ- 
ation happens in the case of KdV and mKdV equations, and more generally, if the 
Hamiltonian is positive definite or if h(u) does not grow too fast at infinity JSjj. 

As a corollary of our theory of the string equation we can prove that all the 
perturbative corrections of the Cauchy problem of KdV are rational functions of 
the solution of the Hopf equation and its derivative. In other words, there exists 
at any order in £ a quasi-Miura transformation mapping the solution of the Hopf 
equation to the asymptotic expansion (for e small) of the solution of KdV. This 
result was proved so far only for solutions of the Hopf equation with never vanishing 
first derivative and for initial data depending on e in a prescribed way, see [4] , . 

Theorem 4. Let u(x,t,e) be the solution of the KdV Cauchy problem 
u t = uu x + e 2 u xxx , u t=Q = (f € P| H s (R), 

s>0 

and let t c > be the critical time of the Hopf equation with same initial data. Then, 
for every t G [0, t c ) and the solution u possesses the asymptotic expansion 

u(x,t,e) = ]TVV(M), (50) 

i>0 

where 

v\x, t) = C(v°(x, t), d x v°(x, t),..., d™*v°(x, t), t), (51) 

for some mi. The functions £ z (i> , v x , . . . ,t),i € N are rational with respect to all the 
variables, and they do not depend on the initial data of the Cauchy problem. 

Proof. From the results of [T3], it follows that for the class of initial data in the 
hypotheses, the solution u admits an asymptotic expansion, for times smaller that 
t c , of the form (fSTTj) . Here, v° is the solution of the Hopf Cauchy problem with same 
initial data, and the smooth functions v l are the unique solutions (in H°°) of the 
semilinear equations 

<•; X (])'•''" '• « l |*=o = 0, i>l. (52) 

j=0 

Given any smooth initial data (p, we know that for any N > 0, the KdV equation 
admits a functional Sf satisfying (|21bl 121c)) . and condition (|21a|) with K = 2N + 2. 
By the results of Section [51 the corresponding string equation is of the form 

x + ut = f'(u) + e 2 'T 3 {f"{u), u x , . . . , «<«>) , (53) 

for some functions Tj, depending rationally on f'(u), . . . , f( k ^(u), and polynomially 
on u x , . . . , u^ 2 ^. As usual, here /' is a local inverse of ip. Note that the Tj are 
independent of the initial data. 

We prove the theorem by induction. The zero-th correction is trivially true by 
taking £°(u°) = v°. Suppose now (|5"Tj) holds for all f for j = 0, . . . , k. 
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Let us restrict our attention to points (x, f) with t G [0, t c ) and such that 
v®(x, i) ^ 0. Indeed, by Theorems [5] and [3] and for smooth initial data, in this 
case all the e-derivatives of (|S"5)) vanish at e = 0. Taking the (2k + 2)-th partial 
derivative of with respect to e, and evaluating at e = 0, we get 

(t - /( fe+2 V)) « fc+1 = R k+1 (f"(v°), . . . , / W (A v° x , . . . , v\ v k . . . ) , (54) 

for some rational function Rk+i- Now since x + v°t = f'(v°), then the quantities 
/W(t;°) can be written as rational functions of t and of the derivatives of v°. Due to 
the induction hypothesis, the same can be said for v l , 1 < i < k and their derivatives. 
Therefore, by (|54j) v k+1 is a rational function of t and of v° and its derivatives. We 
call this function 

It remains to prove that v k+1 (x,t) = £ k+1 (v°(x, t), v x (x, t), . . . , t) also where 
v Q x = 0, provided (x,t) £ R x [0,t c [. This is easily accomplished by noticing that 
f*+i( w o (ar, (a^t),. ..,t) satisfies (|S2J -with i = fc + 1- for any (x,t) G K x [0,t c [. 

In fact ^ fc+1 is analytic at v x = for any value of the jets variables v°,v xx , . . . 
and provided t G [0,t c ), because due to Theorem [3] v k+1 (x, t) is smooth for any 
initial data and for (x,t) elx [0, t c ). Moreover, using the identity u° = v°v x and 
the induction hypothesis on = i,...,k the k + 1-th transport equation (|5^|) 
reduces to an algebraic identity for £ fc+1 which is then satisfied identically on the 
domain of analyticity of £ k+1 . □ 

Conclusions 

We have introduced a deformation of the method of characteristics valid for Hamil- 
tonian perturbations of a scalar conservation law in the small e (dispersionless, 
semiclassical) limit. Our method of analysis is based on the 'variational string equa- 
tion', a functional-differential relation originally introduced by Dubrovin [3] in a 
particular case, of which we have laid the mathematical foundation. Starting from 
first principles, we have constructed the string equation explicitly up to the fourth 
order in perturbation theory and we have shown that the solution of (the Cauchy 
problem of) the Hamiltonian PDE satisfies the appropriate string equation in the 
small e limit and for times smaller than the critical time. Moreover, we applied 
our construction to compute explicitly the first two perturbative corrections of the 
solution of the general Hamiltonian PDE, and we proved for KdV the existence of 
a global quasi-triviality transformation at any order in e. 

Yet many questions remain open about the string equation and its applicability, 
and more in general, about the behaviour of solutions of Hamiltonian PDEs in the 
small e regimes. We plan to address them in future publications. A prominent 
open problem is the rigorous treatment of the dispersive shock phenomenon, the 
behaviour of solutions of Hamiltonian PDEs in the small e regime and close to the 
critical time. According to the Dubrovin universality conjecture, the solution of the 
generic Hamiltonian PDE undergoes a phase transition from a regular behaviour 
to a highly oscillatory behaviour and this transition can be described in term of a 
special solution of the integrable ODE known as PI(2), which is the second member 
of the Painleve I hierarchy. In fact, a formal multiscale expansion centered in a 
neighborhood of the critical time transforms any Hamiltonian PDE to KdV and the 
same multiscale expansion reduces the fourth order string equation of KdV to PI(2). 
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Remarkably, all higher order terms of the string equation do not contribute to the 
expansion, so that the PI(2) correction is beyond all order in perturbation theory 
while based on a fourth order computation. 

Our work provides another argument in favor of the universality conjecture. 
Indeed, the same formal computation shows that the fourth order string equation of 
any Hamiltonian PDE (which we constructed) reduces to PI(2) in that multiscale 
expansion. To make such heuristic considerations into a proof we plan to address 
the validity of the string equation in a neighborhood of the critical time. 

Another problem we plan to address is more algebraic in nature: as we have 
shown, if we consider only functionals with densities polynomial in the jets, then 
not all the Hamiltonian PDEs admit a string equation of order four. Indeed, a non- 
trivial obstruction arises at the fourth order and, presumably, higher order terms 
would impose more constraints on the Hamiltonian. We will investigate whether 
the obstruction(s) can be lifted by considering more general functionals, for example 
rational in the jet variables. 
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